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Abstract 


The  first  and  second  order  temperature  coefficients  of  the  material  constants  c® ,  cf2 ,  4 , 
4,  and  el5  of  dilithium  tetraborate  (U2B4O7)  have  been  measured  over  a  temperature  range  of 

20°C  to  150°C.  An  improved  resonator  method  was  used  to  measure  the  fundamental  zero  mass 
loading  antiresonant  frequencies  of  selected  pure-mode  orientations  of  U2B4O7.  Material 
constants  extraction  was  performed  using  a  linear  least  squares  matrix  method.  The  resulting 
material  constant  curves  were  fit  with  a  third  order  power  series  to  obtain  their  corresponding 
temperature  coefficients.  The  calculated  temperature  coefficients  of  the  material  constants  cft , 
c,2 ,  4  *  4  *  el5  of  U2B4O7  were  used  to  predict  the  zero  mass  loading  antiresonant 
frequencies  of  the  crystal  samples  with  average  errors  of  1.20%  for  thickness  excitation  and 
0.285%  for  lateral  excitation. 
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Determination  of  the  Temperature  Coefficients  of  Selected 
Material  Constants  of  Dilithium  Tetraborate 

Chapter  1  -  Introduction 


This  chapter  provides  the  background  material  necessary  to  understand  the  importance  of 
piezoelectric  crystals  to  present  day  electronic  systems.  Various  applications  of  piezoelectric 
crystals  are  discussed.  The  two  main  categories  of  acoustic  piezoelectric  devices  (bulk  and 
surface  acoustic  wave)  are  presented  and  described.  Problem  Statement  and  Scope  sections  cover 
the  goals  of  this  thesis  research  as  well  as  the  experimental  methods  utilized.  Finally,  the 
Approach  and  Presentation  section  provides  an  outline  for  the  remaining  chapters  of  the  thesis. 

1.1  Background 

Piezoelectric  materials  have  become  essential  to  the  operation  of  many  electronic  systems 
both  in  commercial  and  military  use  due  to  the  fact  that  they  are  natural  transducers.  Their  use  as 
transducers  -  converting  electrical  to  mechanical  energy  or  vice  versa  -  is  an  electric  tap  into  the 
very  desirable  mechanical  properties  of  crystals.  Exploiting  these  desirable  properties  of  crystals, 
such  as  very  high  Q  values,  allows  the  production  of  very  accurate  resonators.  Highly  accurate 
crystal  resonators  are  needed  for  use  in  clocks  and  oscillators  whose  accuracies  are  critical  to  the 
performance  of  many  types  of  electronic  systems  in  the  fields  of  navigation,  communication, 
digital  systems,  space  tracking,  guidance  systems,  radar,  and  electronic  warfare. 

One  of  the  desirable  properties  into  which  piezoelectricity  taps  is  the  high  elastic  constants 
associated  with  crystals.  The  interatomic  bonds  in  a  crystal  hold  the  atoms  rigidly  in  an  orderly 
lattice  structure.  The  extreme  stiffness  of  these  bonds  is  what  gives  crystals  their  characteristic 
hardness  and  correspondingly  their  high  elastic  constant  values.  As  the  value  of  the  elastic 
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constant  increases,  so  does  the  ability  of  the  atoms  in  the  crystal  to  return  to  their  equilibrium 
positions  after  being  perturbed.  The  ability  of  a  material  to  transmit  high  frequency  mechanical 
signals  is  directly  related  to  how  quickly  the  signal  carriers  can  oscillate;  therefore,  the  high  elastic 
constants  associated  with  crystals  allow  them  to  be  used  in  high  frequency  applications  such  as 
filters  and  delay  lines  (Figure  1-1). 


metal  counter 

ELECTRODE 


PIEZOELECTRIC 
MATERIAL 


METAL  TOP 
ELECTRODE 


Fig  1-1.  Example  of  a  high  frequency  delay  line  [1].  The  delay  line  rod  is  not  a 
piezoelectric  crystal.  Only  the  ends  of  the  rod  are  piezoelectric  (transducers). 

This  allows  the  delay  line  rod  to  be  a  material  possessing  more  desirable 
mechanical  properties  than  the  piezoelectric. 

The  “natural  transducer”  characteristic  of  piezoelectric  materials  lends  them  for  use  as  tactile 
sensors  on  a  robotic  finger  or  in  bomb  fuses  to  indicate  when  the  bomb  has  hit  the  target. 
Another  transducer-type  application  is  piezoelectric  transmitters  and  receivers  that  send  out  and 
receive  vibrational  signals.  Piezoelectric  transmitters/receivers  find  use  as  underwater 
communication  systems  or  as  the  ultrasound  machines  used  in  medicine  (Figure  1-2).  In  addition 
to  their  high  Q  values,  piezoelectric  materials  have  other  benefits  such  as:  stability  over  a  wide 
range  of  parameters;  low  noise;  low  power;  small  size;  fast  warm-up;  and  low  life-cost  cycle.  The 
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importance  of  piezoelectric  materials  to  modem  electrical  systems  cannot  be  overemphasized.  As 
our  knowledge  about  piezoelectric  materials  grows,  so  does  our  ability  to  design  and  build  more 
advanced  systems. 


Fig  1-2.  Example  of  a  piezoelectric  transmitter/receiver  [1].  The  matching  layer 
serves  to  match  the  acoustic  impedances  of  the  piezoelectric  crystal  and  the 
material  in  which  it  transmits. 


There  are  two  general  classifications  of  acoustic  piezoelectric  devices:  bulk  acoustic  wave 
(BAW)  and  surface  acoustic  wave  (SAW),  as  shown  in  Figure  1-3.  The  difference  between  the 
two  types  of  devices  lies  in  the  type  of  elastic  waves  generated  in  the  crystal.  The  BAW  device 
has  “sending”  electrodes  that  create  acoustic  waves  that  travel  through  the  bulk  of  the  material 
and  are  picked  up  by  “receiving”  electrodes  that  can  be  on  either  side  of  the  crystal.  Types  of 
BAW  devices  include  resonators,  filters,  transducers  from  audio  frequencies  to  the  UHF  band, 
and  oscillators  [2].  The  SAW  devices  have  “sending”  electrodes  that  create  acoustic  waves  which 
travel  on  the  surface  of  the  crystal  and  are  picked  up  on  the  same  surface  by  the  “receiving” 
electrodes.  Types  of  SAW  devices  include  resonators,  delay  lines,  filters,  and  correlators  [2]. 
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BAW  Device  SAW  Device 

Fig  1-3.  Examples  of  a  bulk  acoustic  wave  (BAW)  and  a  surface  acoustic  wave 

(SAW)  devices  [3]. 

1.2  Problem  Statement 

With  the  importance  of  piezoelectric  materials  fully  realized,  it  is  understood  that 
knowledge  about  the  material  properties  of  piezoelectric  materials  is  necessary  in  order  to  develop 
and  advance  the  electrical  systems  that  use  them.  Quartz  crystals  have  been  the  most  studied  and 
most  used  piezoelectric  crystals  to  date;  however,  studies  of  dilithium  tetraborate  (U2B4O7)  have 
demonstrated  that  this  material  might  do  a  better  job  than  quartz  in  many  applications  [2].  Given 
the  correct  temperature  coefficients  of  the  relevant  material  constants  (elastic  constant, 
piezoelectric  constant,  dielectric  constant),  the  resonator  frequencies  of  a  material  may  be 
predicted  at  any  temperature.  The  temperature  coefficients  of  the  material  constants  of  U2B4O7 
published  to  date  have  not  been  able  to  accurately  predict  the  temperature  behavior  of  its 
resonator  frequencies  [4].  In  order  to  fully  realize  the  abilities  of  dilithium  tetraborate,  accurate 
knowledge  of  the  temperature  coefficients  of  its  material  constants  must  be  obtained.  The  goal  of 


1-4 


this  thesis  research  is  to  determine  the  temperature  dependence  of  selected  material  constants  of 
dilithium  tetraborate. 

1.3  Scope 

Samples  of  dilithium  tetraborate  were  obtained  from  Dr.  John  A.  Kosinski,  Army  Research 
Laboratory,  Fort  Monmouth,  NJ.  It  was  necessary  to  design  and  build  a  test  fixture  to  hold  the 
1^2^407  crystal  samples.  The  resulting  test  fixture  serves  the  dual  purpose  of  holding  the 
crystals  and  stabilizing  the  temperature  of  the  crystals.  The  antiresonant  frequencies  of  the 
samples  were  measured  over  a  range  of  temperatures.  The  resulting  frequency  dependences  of 
temperature  were  modeled  by  a  third  order  power  fit.  Using  previously  published  values  for  the 
zero  mass-loading  room  temperature  material  constants  of  U2B4O7  [2],  the  mass-loading  effect 
was  removed  from  the  fitted  frequency  versus  temperature  curves.  The  resulting  set  of  frequency 
versus  temperature  curves  were  used  in  a  linear  least  squares  extraction  of  the  selected  material 
constants  of  U2B4O7  over  a  range  of  temperatures.  Temperature  coefficient  curves  for  the 
selected  material  constants  of  U2B4O7  were  then  constructed  from  the  results  of  the  material 
constants  extraction.  Results  are  compared  with  other  results  published  in  the  literature. 
Suggestions  for  further  work  are  also  presented. 

1.4  Approach  and  Presentation 

The  goal  of  this  research  was  to  obtain  accurate  values  for  the  temperature  coefficients  of 
selected  material  constants  of  dilithium  tetraborate.  Chapter  2  presents  basic  theory  beginning 
with  the  conceptual  fundamentals  of  piezoelectricity  along  with  a  review  of  tensor  notation  in  the 
context  of  the  piezoelectric  constitutive  equations.  Chapter  2  concludes  with  a  description  of  the 
“resonator  method”  used  in  this  thesis  as  well  as  a  definition  of  the  temperature  coefficients  whose 


1-5 


determination  is  the  ultimate  goal  of  this  research.  Chapter  3  presents  a  detailed  derivation  of  the 
eigenvalue  expressions  used  to  obtain  pure-mode  solutions  to  the  piezoelectric  wave  equation. 
Chapter  4  describes  the  experimental  setup  as  well  as  the  methods  used  to  extract  the  temperature 
coefficients  of  the  material  constants.  Chapter  5  contains  the  results  of  this  research  as  well  as  a 
comparison  with  other  published  results.  Chapter  6  presents  conclusions  drawn  from  this 
research  and  suggestions  for  further  investigation. 
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Chapter  2  -  Fundamentals 


This  chapter  provides  an  overview  of  the  theory  of  piezoelectricity.  Basic  concepts  of 
piezoelectricity  are  presented  as  well  as  the  physical  models  and  mathematical  formulas  on  which 
piezoelectric  phenomena  are  based.  A  review  of  tensor  notation  is  presented  in  the  context  of  the 
piezoelectric  constitutive  relations.  The  experimental  method  used  in  this  thesis  for  determining 
the  material  constants  of  a  piezoelectric  crystal,  the  “resonator  method,”  is  reviewed.  Finally, 
there  is  a  brief  discussion  of  the  temperature  dependence  of  material  constants. 

2.1  Concepts 

Piezoelectric  properties  are  found  in  many  crystals.  A  piezoelectric  material  will  deform 
when  an  electric  field  is  applied  across  it,  and  conversely,  an  electric  field  will  be  generated  across 
a  piezoelectric  material  when  it  is  stressed.  The  cause  of  piezoelectricity  is  rooted  in  the 
symmetry  of  a  material.  Materials  lacking  a  center  of  symmetry  (non-centrosymmetric)  are 
piezoelectric.  Figure  2-1  shows  a  material  that  possesses  trifold  symmetry.  The  electric  dipole 
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Fig  2-1.  Crystal  possessing  trifold  symmetry  [1]. 
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moments  (arrows)  add  up  to  zero  in  the  unstressed  case;  however,  in  the  stressed  case  the  dipole 
moments  add  up  to  give  a  net  electric  dipole  moment  upwards.  The  crystal  structure  of  U2B4O7 
is  shown  in  Figure  2-2. 


Fig  2-2.  The  structure  of  dilithium  tetraborate  as  viewed  along  the  x3-axis  [5].  The 
black  spheres  represent  boron  atoms,  the  white  spheres  represent  oxygen  atoms,  and 
the  shaded  spheres  represent  lithium  atoms.  The  square  represents  a  unit  cell. 

The  basic  structural  unit  of  dilithium  tetraborate  consists  of  two  non-planar  six  membered 
rings  combined  in  one  group  [5].  Because  the  group  is  slightly  twisted,  it  has  no  true  mirror- 
plane  symmetry  in  the  x3  direction;  thus,  dilithium  tetraborate  is  a  piezoelectric  material.  There  is, 
however,  a  four-fold  symmetry  about  the  x3-axis,  which  places  L^B 4O7  in  the  4mm  crystal  class. 
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A  crystal  can  be  modeled  as  a  damped  mass-spring  system  in  which  the  atoms  (mass)  are 
held  together  in  a  lattice  by  chemical  bonds  (springs)  which  have  some  form  of  damping 
associated  with  them  (dashpot).  A  mass-spring  system  behaves  as  an  oscillator,  so  one  can  offer 
another  equivalent  model  for  the  crystal.  This  model  is  an  electrical  resonator,  or  RCL  circuit. 
Both  models  are  shown  in  Figure  2-3. 


Fig  2-3.  Mechanical  model  of  crystal  and  equivalent  electrical  model  [3]. 

For  comparison,  the  mass  in  the  mechanical  model  corresponds  to  the  inductor  in  the 
electrical  model.  Both  the  mass  and  the  inductor  are  a  resistance  to  change  (mass  resists  change 
in  motion,  inductor  resists  change  in  current).  The  spring  in  the  mechanical  system  corresponds 
to  the  capacitor  in  the  electrical  system.  Both  the  spring  and  the  capacitor  store  energy  (the 
spring  stores  mechanical  energy,  the  capacitor  stores  electrical  energy).  The  dashpot  corresponds 
to  the  resistor.  Both  are  resistant  to  “motion”,  whether  it's  mechanical  motion  or  the  electrical 
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“motion”  of  current,  and  both  represent  loss.  Solving  the  differential  equations  governing  these 
systems,  one  would  obtain  solutions  in  the  form  of  waves.  In  Chapter  3  the  governing  differential 
equations  are  solved  to  yield  wave  solutions  that  are  dependent  on  the  material  constants  as  well 
as  the  orientation  of  the  sample  relative  to  the  crystallographic  axes.  The  waves  in  a  piezoelectric 
material  can  be  either  transverse  (shear),  longitudinal,  or  a  combination  of  both. 

From  the  RCL  model  of  a  piezoelectric  resonator  it  is  possible  to  obtain  the  equivalent 
circuit  for  the  bulk  acoustic  wave  resonator  shown  in  Figure  1-3.  The  equivalent  circuit  is  shown 
in  Figure  2-4.  The  resistor  (R),  capacitor  (C),  and  inductor  (L)  represent  electrical  equivalent 
model  for  the  vibrating  crystal  and  the  capacitance  Cq  is  the  capacitance  resulting  from  the 
attachment  of  electrodes  to  the  crystal  in  order  to  create  the  acoustic  waves.  Since  the  electrical 
components  of  the  equivalent  circuit  correspond  to  the  mechanical  components  of  the  crystal,  the 
characteristic  resonances  can  be  used  to  measure  the  material  constants  of  piezoelectric  crystal. 


Co 


Fig  2-4.  Equivalent  circuit  for  a  BAW  resonator  [3].  Where  C  is  the  capacitance, 
L  is  the  inductance,  and  R  is  the  resistance  associated  with  the  electrical  equivalent 
model  of  the  crystal.  Co  is  capacitance  added  as  a  result  of  the  deposition  of 
electrodes  onto  the  crystal  surface. 
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The  electric  and  mechanical  properties  of  a  crystal  are  coupled  together  through  the 


piezoelectric  constitutive  relations  [1]: 


T=c?  S- - e ..  E. 

>  '-y  uj  ^ j 
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(2.1a) 

(2.1b) 


where:  T  is  stress;  cE  is  elastic  stiffness  at  a  constant  electric  field;  S  is  strain;  e  is  piezoelectric 
constant;  E  is  electric  field;  D  is  electric  displacement;  £s  is  dielectric  constant  at  a  constant  strain. 
Stress  (T)  is  the  force  per  unit  area  that  is  applied  to  a  solid.  Strain  (S)  is  the  fractional  particle 
displacement  that  occurs  due  to  stress.  The  elastic  constant  (cE)  is  a  measure  of  how  strongly  the 
atoms  of  a  material  are  bound  to  each  other  at  a  constant  electric  field  value.  The  piezoelectric 
constant  (e)  is  a  measure  of  how  strongly  the  electrical  and  mechanical  properties  of  a  solid  are 
coupled.  The  elastic,  piezoelectric,  and  dielectric  constants  are  the  material  properties  that  govern 
the  piezoelectric  effects  of  a  material.  These  material  constants  are  tensors,  and  the  piezoelectric 
effects  seen  in  crystals  are  generally  anisotropic.  In  the  design  of  piezoelectric  devices  it  is 
necessary  to  know  these  material  constants  to  a  high  degree  of  accuracy. 


2.2  Tensor  Notation 

The  piezoelectric  constitutive  relations  may  be  written  in  a  matrix  form  [3]: 
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The  tensor  values  represent  the  fact  that  the  material  is  anisotropic,  but  they  also  represent  the 
fact  that  there  are  two  ways  to  apply  stress  to  a  material.  Stress  can  be  applied  either 
longitudinally  (normal)  or  transversely  (shear)  to  a  material  surface.  The  components  of  the  stress 
tensor  are  shown  pictorially  in  Figure  2-5.  The  first  subscript  denotes  the  coordinate  axis  normal 


to  a  stressed  plane,  and  the  second  subscript  denotes  the  axis  along  which  the  stress  is  applied. 
The  diagonal  element  of  the  stress  tensor,  for  example  Txx,  represents  a  stress  in  the  x  direction 
due  to  a  force  in  the  x  direction.  The  off  diagonal  element  of  stress,  for  example  Txy,  represents  a 
stress  on  the  x  plane  due  to  a  force  in  the  y  direction.  Equation  (2.2)  makes  use  of  Voigt’s 
prescription  which  takes  into  account  the  inherent  symmetry  of  crystals.  The  use  of  Voigt’s 
prescription  reduces  the  rank  of  a  symmetric  tensor  by  one  half,  thereby  greatly  simplifying  the 
problem.  Voigt’s  prescription  is  as  follows  [6]: 

ij  or  kl  11  22  33  23,32  31,13  12,21 

X  or  p.  1  2  3  4  5  6 
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where  the  subscripts  ij  or  kl  may  be  replaced  with  the  corresponding  value  of  X  or  (X.  The  values 
Ti,  T2,  and  T3  correspond  to  T^,  Tyy,  and  Tzz  which  are  the  longitudinal  stresses  in  the  x,  y,  and  z 
directions  respectively^  the  values  T4,  T5,  and  Tg  correspond  TyZ  or  TZy,  Tzx  or  T x/,  TXy  or  TyX 
which  are  shear  stresses  about  the  x,  y,  and  z,  axis  respectively.  Equation  (2.2)  gives  nine 
coupled  equations  that  relate  stresses,  strains,  piezoelectric  constants,  elastic  constants,  electric 
fields,  dielectric  constants,  and  electric  displacements  in  different  directions. 

It  appears  necessary  to  obtain  values  for  eighty  one  material  constants  in  order  to 
completely  characterize  a  piezoelectric  material!  In  practice,  however,  the  number  of  independent 
entries  in  the  material  constants  matrix  in  Equation  (2.2)  is  greatly  reduced  due  to  symmetries 
arising  from  the  crystal  structure.  Dilithium  tetraborate,  for  example,  belongs  to  symmetry  class 
4mm  crystal  with  two  independent  values  for  dielectric  constant,  three  independent  values  for 
piezoelectric  constant,  and  six  independent  values  for  elastic  constant  [2].  The  dielectric  constant 
matrix  of  dilithium  tetraborate  is  as  follows  [7]: 

£n  0  0 

[e]=  0  £,,  0  (2.3a) 

0  0  £33 

The  piezoelectric  constant  matrix  of  dilithium  tetraborate  is  as  follows  [7]: 

‘  0  0  0  0  e15  O' 

[e]=  0  0  0  e15  0  0  (2.3b) 

e31  e31  «33  0  0  0 

The  elastic  constant  matrix  of  dilithium  tetraborate  is  as  follows  [7]: 
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[c]= 


(2.3c) 


cu  Cj2  cI3  0  0  0 

C12  C11  C13  0  0  0 

c13  c13  c33  0  0  0 

0  0  0  C44  0  0 

0  0  0  0  044  0 

0  0  0  0  0  % 

So  the  characterization  of  U2B4O7  is  reduced  to  a  problem  of  finding  eleven  independent 
variables. 

2.3  Measurement  Procedure  (The  “Resonator  Method”) 

A  measurement  technique  presented  by  Kosinski  [2]  has  advantages  over  other  methods 
arising  from  the  exploitation  of  pure-mode  excitation  of  the  material  combined  with  an 
engineering  methodology  for  pure-mode  calculations.  The  measurement  technique  is  described  in 
detail  in  Chapter  4;  however,  a  basic  overview  is  presented  here. 

The  technique  employs  the  “resonator  method”  in  which  a  sample’s  characteristic 
antiresonance  frequencies  are  measured.  Multiple  odd  harmonics  of  the  fundamental 
antiresonance  frequencies  are  measured  in  order  to  increase  the  accuracy  of  the  results.  The 
characteristic  antiresonance  frequencies  depend  on  the  material  constants  and  the  geometry  of  the 
device.  The  crystals  used  in  the  technique  are  cut  into  discs  that  measure  0.2  mm  in  thickness  and 
14  mm  in  diameter  [2].  Consider  the  experimental  setup  in  Figure  2-6.  This  setup  can  be  used  to 
produce  longitudinal  waves  (pure-mode  device),  shear  waves  (pure-mode  device),  or  a 
combination  of  longitudinal  and  shear  waves  (non  pure-mode  device). 

After  determining  the  antiresonance  frequencies  of  the  crystal  plates,  the  stiffness 
eigenvalues  ( c )  may  be  determined  using  the  following  equation  [2]: 
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/„= 


M 


r-i/2 


2(2  h) 


(2.4) 


where  fA  is  the  antiresonance  frequency,  M  is  the  mode  of  vibration,  2 h  is  the  thickness  of  the 
plate,  and  p  is  the  mass  density  of  the  crystal.  The  stiffness  eigenvalues  are  related  to  the  stiffened 
elastic  constants  by  the  following  differential  equation  (wave  equation): 

_  a2- 


^dx2 


+  cu  =  0 


(2.5) 


where  u  is  particle  displacement,  and  is  the  stiffened  elastic  constant.  Equation  (2.5)  and  the 


eigenvalue  expressions  are  derived  in  Chapter  3. 


T 

2h 


Fig  2-6.  Experimental  setup.  BAW  resonator  with  thickness  2h. 

The  eigenvalue  expressions  for  the  stiffness  eigenvalues  of  a  piezoelectric  plate  form  an 
overdetermined  set  of  equations  relating  the  unknown  material  constants  to  the  known  stiffened 
elastic  constants.  The  antiresonant  frequency  data  from  different  crystal  orientations  as  well  as 
different  excitation  methods  is  combined  with  the  eigenvalue  expressions  in  a  Unear  least  squares 


2-9 


extraction  process.  The  result  of  the  linear  least  squares  extraction  is  the  desired  material 
constants. 

2.4  Temperature  Dependence 

It  is  desirable  to  produce  devices  that  give  the  same  performance  no  matter  what  the 
operating  temperature.  Therefore,  the  temperature  dependence  of  the  material  constants  is  an 
important  characteristic  to  investigate.  The  temperature  dependence  of  a  function  is  modeled  by  a 
power  series  expansion  about  a  reference  temperature.  Using  a  reference  temperature  of  25°C, 
the  resonant  frequency  as  a  function  of  temperature  is  represented  in  the  following  equation: 


/(0)  =  /(0=25°C> 


l  +  7;(0-25)+r2(0-25)2+  ••• 
-+Tn(Q-25j 


(2.6) 


The  value  0  is  the  temperature  and  the  values  Ti,  T2,...,  Tn  are  the  first  through  nth  order 
temperature  coefficients.  In  the  case  of  some  important  piezoelectric  materials  which  have  zero 
temperature  dependence  of  frequency  or  delay  along  specific  cuts,  the  temperature  coefficients  in 
the  power  series  expansion  are  very  small.  As  the  temperature  coefficients  get  very  small,  the 
power  series  expansion  of  frequency  or  delay  loses  its  dependency  on  temperature. 

Since  the  resonator  method  uses  measurements  of  antiresonant  frequencies  to  extract  the 
material  constants,  when  the  antiresonant  frequencies  change  with  temperature  so  will  the  material 
constants.  The  plate  thickness  (2 h),  mass  density  (p),  and  dielectric  constant  (e)  are  also  functions 
of  temperature.  Because  they  are  functions  of  temperature,  2 h,  p,  and  e  also  have  their  own 
power  series  expansions.  If  the  correct  temperature  coefficients  for  each  of  the  power  series 
expansions  were  know,  it  would  be  possible  to  predict  the  values  of  the  material  constants  at 
every  temperature.  In  order  to  investigate  the  temperature  dependence  of  the  material  constants, 
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it  is  necessary  to  determine  the  temperature  dependence  of  frequency,  plate  thickness,  mass 
density,  and  dielectric  constant.  Fortunately,  the  temperature  dependences  of  plate  thickness  and 
mass  density  are  the  result  of  simple  thermal  expansion.  The  coefficients  of  thermal  expansion 
have  been  determined  by  other  researchers  [8-14]. 

It  is  desired  to  have  materials  that  possess  little  or  no  temperature  dependence  for  certain 
device  characteristics  such  as  frequency  or  delay.  With  accurate  knowledge  of  the  room 
temperature  values  for  the  material  constants  and  their  temperature  coefficients  it  will  be  possible 
to  predict  the  material’s  behavior  at  any  orientation  and  at  any  temperature,  thus  allowing  full 
exploitation  of  the  material. 
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Chapter  3  -  Theory 

This  chapter  provides  the  reader  with  a  detailed  derivation  of  acoustic  plane  wave 
propagation  in  piezoelectric  materials.  Once  a  solution  for  acoustic  wave  propagation  is  obtained, 
generalized  rotation  matrices  are  presented  so  that  wave  propagation  in  any  direction  through  a 
material  may  be  examined.  Applying  conditions  for  pure-mode  solutions  allows  the  determination 
of  a  set  of  crystal  orientations  that  propagate  pure-mode  waves.  Selected  crystal  orientations  that 
propagate  pure-mode  waves  are  used  in  the  experimental  procedure  described  in  Chapter  4. 

3.1  Acoustic  Wave  Propagation  in  Piezoelectric  Materials 

In  order  to  determine  the  crystal  orientations  that  will  yield  pure-mode  waves  when 
excited,  it  is  necessary  to  examine  how  acoustic  waves  travel  in  a  piezoelectric  crystal.  The 
derivation  in  this  section  is  taken  from  Tiersten’s  paper  on  thickness  vibrations  in  piezoelectric 
plates  [15].  Some  steps  were  filled  in  where  appropriate. 

Consider  an  unrotated  flat  piezoelectric  plate  of  infinite  extent  (Figure  3-1).  Using  the 
rectangular  Cartesian  coordinate  system,  let  the  piezoelectric  plate  be  infinite  in  the  xi  and  x3 

directions,  and  choose  the  faces  of  the  plate  to  be  defined  at  X2  =  ±h. 

♦X2 


Fig  3-1.  Infinite  flat  piezoelectric  plate  [2], 
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Differentiation  with  respect  to  a  spatial  coordinant  will  be  denoted  by  preceding  the  differentiation 
index  by  a  comma.  The  dot  notation  for  differentiation  with  respect  to  time  and  the  summation 


convention  for  repeated  indices  are  also  employed. 

The  system  of  governing  equations  for  this  infinite  piezoelectric  plate  includes  the  stress 
equations  of  motion: 

TiU=puj  (3.1) 


The  charge  equation  of  electrostatics: 


D  =  0 

M 


The  mechanical  strain  displacement  relations: 

Ski  =l(Uk,l  +Ul,k) 


(3.2) 


(3.3) 


The  equation  of  electrostatics: 


Ek  -“9,1 


(3.4) 


The  linear,  piezoelectric  constitutive  relations: 


Tij  —  C ijkl  '  Skl  ekij 

■Ek 

(3.5a) 

A  =eit  *  A  +  em 

■Ski 

(3.5b) 

where  all  symbols  in  equations  3.1  through  3.5  have  been  defined  in  Chapter  2  except  for  (p, 
which  is  electric  potential. 

An  alternating  electric  potential  is  applied  to  the  faces  of  the  plate.  Assuming  no  external 
forces  are  applied  to  the  plate  face,  the  boundary  conditions  for  the  infinite  plate  are: 


II 

o 

at 

*2  =±H 

(3.6) 

<P=±<Po*"i<°' 

at 

x2  =±h 

(3.7) 
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where  cpo  is  a  constant,  and  ©is  angular  frequency. 

Now  Equations  (3.3)  and  (3.4)  may  be  substituted  into  Equation  (3.5a)  to  yield: 

Tij  =CijklUk,l+ekij(?,k  (3.8a) 

Equations  (3.3)  and  (3.4)  may  be  substituted  into  Equation  (3.5b)  to  yield: 

A  ~eiklUk,l  ~  (3.8b) 


Since  there  are  no  boundaries  in  the  x,  and  x3  directions  and  the  applied  voltage  is  independent  of 
the  xi  and  x3  directions,  all  derivatives  with  respect  to  these  directions  will  vanish.  Substituting 
Equation  (3.8a)  into  Equation  (3.1)  gives: 


C2  jk2Uk,22  +  e22j(?  ,22  ~ 

Substituting  Equation  (3.8b)  into  Equation  (3.2)  yields: 

e2k2Uk,22  ~  ^22*P,22  =  ® 

Now  substituting  Equation  (3.8a)  into  Equation  (3.6)  gives  the  boundary  condition: 

C2jk2Uk,2  ~^e22jty,2  ~  ®  ^ 


x2  =±h 


along  with  the  other  boundary  condition  (3.7).  Rearranging  Equation  (3.9b)  gives: 


<P,22  = 


'2k2 


’.22  —  s  Uk,22 
E22 


Substituting  Equation  (3.11)  into  Equation  (3.9a)  yields: 


(3.9a) 


(3.9b) 


(3.10) 


(3.11) 


C2jk2Uk,22  +  e22  j 


—  U 

S  uk.  22 

\E22  J 


-  p  Uj  =  0 


(3.12) 


Assuming  an  e'1  time  dependence  and  combining  terms.  Equation  (3.12)  becomes: 

C2jk2Uk,22  +  P®  Uj  =  ® 


(3.13a) 
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where: 


(3.13b) 


—  _  E  .  e22je22k 

C2jk2  ~  C2jk2  ■*  ~S 
^22 

The  value  c2jk2  is  the  stiffened  elastic  constant.  Equation  (3.13)  is  the  governing  differential 

equation  for  thickness  vibrations  in  an  infinite,  flat  piezoelectric  plate.  Solutions  to  this  wave 
equation  are  of  the  form: 


Uj  =  Aj  sin(rp:2) 

The  value  r|  is  the  wave  number  defined  by  Equation  (3.15). 


(3.14) 


(3.15) 


The  value  X  is  the  acoustic  wavelength.  Substituting  Equation  (3.14)  into  the  governing 
differential  equation  (3.13)  gives: 


C2jk2 


pco 


*1 


2 

2 Sjk 


\Ak=0 


where  5jk  is  the  Kronecker  delta.  For  non-trivial  solutions: 


C2jk2  jk  ~  ® 


where: 


_  P®' 


(3.16) 


(3.17a) 

(3.17b) 


Equation  (3.17)  is  the  governing  eigenvalue  equation  for  acoustic  waves  in  a  piezoelectric  plate. 
It  is  a  cubic  equation  in  c  that  yields  three  positive,  real  roots  (c(1),c(2),c(3))  called  the  stiffness 

eigenvalues.  The  eigenvectors,  Ak,  may  then  be  found  using  Equation  (3.16). 

For  the  infinite  flat  piezoelectric  plate  of  Figure  3-1,  acoustic  resonance  occurs  when  the 
thickness  of  the  plate  equals  an  odd  integer  multiple  of  acoustic  half-wavelengths.  This  condition 
is  given  in  Equation  (3.18). 
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where  M  =  1,3,5... 


(3.18) 


When  this  condition  is  met,  the  acoustic  plane  waves  traveling  in  the  x2  direction  will  be 
completely  reflected  at  the  plate  face.  When  the  crystal  resonator  of  Figure  2-6  is  connected  to  a 
transmission  line  circuit,  this  complete  reflectance  corresponds  to  an  open-circuit  (infinite 
impedance)  condition.  The  frequencies  at  which  the  crystal  resonator  acts  like  an  open  circuit  are 
the  antiresonant  frequencies  (fA0)  of  the  circuit.  In  order  to  determine  the/A0’s,  Equation  (3.18)  is 
solved  for  X  and  substituted  into  Equation  (3.15)  to  yield: 


(3.19) 


Equation  (3.19)  is  then  substituted  into  Equation  (3.17b)  to  give: 


c  =  p 


(3.20) 


For  the  fundamental  zero  mass  loading  antiresonant  frequency  (/Ao),  Equation  (3.20)  reduces  to 
the  following: 

c  =  f>(2(2h)Cf  (3.21) 

Equation  (3.21)  is  used  to  obtain  the  stiffness  eigenvalues  from  measurements  of  p,  2h,  and 
using  unelectroded  crystal  (zero  mass  loading). 

3.2  Coordinate  Transformations  and  Symmetry  Conditions 

The  remainder  of  Chapter  3  is  based  on  the  method  presented  by  Kosinski  to  determine 
crystal  orientations  that  propagate  pure- mode  waves  [2].  Once  pure-mode  orientations  are 
determined,  it  is  possible  to  construct  a  sample  set  of  crystal  orientations.  Use  of  the  sample  set 
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with  the  resonator  method  utilizes  the  advantages  of  pure-mode  waves  in  the  determination  of  the 
material  constants. 


Symmetry  considerations  allow  the  governing  eigenvalue  equation  (3.17a)  to  be  written 


using  Voigt’s  prescription  described  in  Chapter  2: 


(^"26  ) 
fe)  -c 

fe) 


(3.22) 


Equation  (3.13b)  may  be  written  with  Voigt’s  prescription  as: 


(3.23) 


In  Equations  (3.22)  and  (3.23)  the  prime  (')  indicates  a  transformation  from  the  crystallographic 
coordinate  system  to  the  coordinate  system  of  the  piezoelectric  plate.  In  matrix  notation,  the 
transformations  are  as  follows  [16]: 


(e5). 

(3.24) 

=WWM' 

(3.25) 

>*)'" 

=  [M][c£][M]r 

(3.26) 

The  matrix  [a]  consists  of  direction  cosines  for  a  doubly  rotated  crystal  plate.  The  matrix  [M]  is  a 
rotation  matrix  formed  from  matrix  [a].  The  superscript  “7"  denotes  the  transpose  of  a  matrix. 

The  crystal  orientations  presented  in  this  thesis  are  based  on  the  standard  IEEE  notation 
(YXwl)<|),0  for  a  doubly  rotated  plate  as  shown  in  Figure  3-2  [17].  The  first  letter  in  the  IEEE 
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*3 

Fig  3-2.  Doubly  rotated  plate  with  (YXwl)<j),0  notation  [17].  The  unprimed  coordinates 
Xi,  X2,  and  X3,  correspond  to  the  crystallographic  axes.  The  double  primed  coordinates 
are  the  doubly  rotated  axes  corresponding  to  the  actual  crystal  cut. 

notation  (YXwl)<|>,0  denotes  the  crystal  axis  which  is  normal  to  the  unrotated  plate  surface  (see 
Figure  3-1).  Thus,  plate  thickness  of  the  unrotated  plate  is  along  the  y-axis  (X2-axis).  The 
second  letter  in  (YXwl)<(),0  denotes  the  direction  of  the  length  of  the  plate;  the  unrotated  plate 
length  is  along  the  x-axis  (Xraxis).  The  width  of  the  crystal  plate,  then,  is  along  the  z-axis  (X3- 
axis).  The  “w”  denotes  the  first  right-handed  rotation  about  the  width  axis  (X3-axis)  of  <|>  degrees. 
This  rotation  results  in  a  new  coordinate  system  for  the  crystal  plate  of  x  i,  x  2,  and  x  3  (x  3  equals 
X3).  The  “1”  denotes  the  second  right-handed  rotation  about  the  singly  rotated  length  axis  (x  r 
axis)  of  0  degrees.  The  result  is  a  doubly  rotated  plate  with  rotated  crystal  axis  of  x  i,  x  2,  and 
x”3. 
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Based  on  the  standard  IEEE  rotation  convention  of  Figure  3-2,  the  direction  cosine  matrix 


is  [16]: 


[a]  = 


cos<()  sin(|)  0 

-sin<|)cos0  cos<()cos0  sin0 
sin<}>sin0  -cos(|)sin0  cos0 


(3.27) 


The  rotation  matrix,  [M],  is  then  given  by  [16]: 


[M]  = 


2 

a\\ 

2 

a\2 

2 

a\3 

2  ana\3 

2anan 

2anan 

Jl 

a2\ 

2 

a22 

2 

a23 

2&22a23 

2&23a2\ 

2a2\a22 

2 

a3 1 
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an 

2 

a33 

2a32&33 

2^33^31 

2a3ld32 

(3.28) 

2ia31 

a22a32 

a23a33 

a22a33  a23a32 

^21^33  +a23a31 

a22a31  +a2\a32 

[3iail 

ayia\2 

a33a\3 

a\2a33  ^13^32 

a\3a3\  +  a\\a33 

a\\a32  +  a\2a3\ 

na2\ 

a\2a22 

a\3a23 

a\2a23  a\3a22 

a\3a2\  ^a\\a23 

a\\a22  J'r(X\2a2\_ 

The  transformation  equations  (3.24)  through  (3.26)  are  applied  to  the  material  constant  matrices 
(2.3a)  through  (2.3c).  The  result  is  a  set  of  analytic  expressions  that  are  used  to  define  the 


elements  of  the  eigenvalue  equation  (3.22)  in  terms  of  untransformed  material  constants  and 


rotation  angles  [2]: 


(£22)  =  cos20(£f,)  +  sin20(£33) 

(3.29) 

/ 

(e22)  =  cos20  sin0(2c15  +  e31)+sin30(e33) 

(3.30) 

/ 

(e24)  =cos30(c15)+cos0sin20(e33  -el5  -e31) 

(3.31) 

(g26 )  =  ® 

(3.32) 

(4)  =cos40[4  - (2 sin 2 (() cos2 <t>)(4  “(4  +46))] 

+2  cos2  0  sin 2  0 (4  +  24 )  +  sin4  0  (4 ) 

(3.33) 
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/ 

(4)  =  cos3 0 sin 0 [(4  +24) -4  +  (2 sin 2 (]) cos2 <() )(4  “(4  +  24)) 


+cos0sin30(4 -(4  +  24))  (3.34) 

/ 

(4)  =  cos30[sin<|)cos<|)(cos2(|)  -sin2<|))(4  ~(4  +  24))  (3.35) 

/ 

(4 )  =  cos4  0(4)  +  cos2  0  sin 2  0  [(4  +  4 ) -  2(4  +  4 ) 

— (2sin2<|)cos2(|))(4  -  (4  +  2c6£6 ))]  +  si n 4  0 (c^ )  (3.36) 

/ 

(4)  =  cos20sin0  sin<])cos<|)(sin2<|) -cos2(|))(4  “(4 +2c6ei))  (3.37) 

(4)  =cos20[4  +  (2sin2<|>cos24>)(4  - (4  +  24 ))]  +  sin 2 0 (4 )  (3.38) 


3.3  Pure-Mode  Solutions  and  Sample  Set  Selection 


Now  that  Equation  (3.22)  may  be  defined  in  terms  of  untransformed  material  constants 
and  rotation  angles,  it  is  possible  to  determine  which  crystal  orientations  will  propagate  pure¬ 
mode  waves.  Pure-mode  solutions  to  Equation  (3.22)  are  obtained  whenever  two  or  more  of  the 


off  diagonal  stiffened  elastic  constant  terms 


/  /  / 
(4)  .  (c26)  ,  (c46) 


are  zero  [2].  When  two  of 


the  off  diagonal  elements  are  zero,  there  is  a  single  distinct  root.  When  all  three  of  the  off 

diagonal  elements  are  zero,  there  are  three  distinct  roots.  In  analyzing  the  zeros  of  the  off 

'  /  ' 

diagonal  elements,  it  is  necessary  to  examine  the  zeros  of  (4  )  ,  (e2X)  ,  and  (e2ll )  . 

Examining  Equations  (3.30)  -  (3.32),  (3.34),  (3.35),  and  (3.37),  it  is  apparent  that  certain 

'  /  ' 

rotation  angles  will  produce  zero  values  for  (4  )  ,  (e2X  )  ,  or  (e2il )  independent  of  the  material 
constant  values.  These  zeros  are  classified  as  global  symmetry  class  zeros  [2].  The  other  class  of 
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zero  loci  occurs  when  combinations  of  the  material  constants  sum  to  produce  a  zero  result  for 

>  ,  > 

(c£)  ,  (e2}  )  ,  or  (e2(1)  •  These  zeros  are  classified  as  material-specific  zeros  [2].  For  dilithium 

/ 

tetraborate,  only  {c24 )  has  material-specific  zero  loci  [2].  The  piezoelectric  contribution  to 

/  / 

(c24)  is  small,  therefore  it  shares  the  same  material-specific  zero  loci  as  (c24 )  [2].  At  <f>  =  0°, 

/ 

material-specific  zeros  for  (c24 )  occur  at  0  =  22.1°,  157.3°,  202.7°,  and  337.3°  [2].  As  (j) 

increases  from  0°  to  45°,  the  zeros  shift  until  at  <|)  =  45°  they  are  located  at  0  =  32.4°,  147.6°, 
212.4°,  and  327.6°  [2], 

Finding  the  pure-mode  solutions  for  dilithium  tetraborate  then  becomes  a  problem  of 

/  /  / 

finding  which  rotation  angles  (<J)  and  0)  produce  zeros  for  at  least  two  of  (c24)  ,  (c26)  ,  and  (c46) 

simultaneously.  The  symmetry  of  class  4mm  materials  limits  the  rotation  angles  to  0°<<t><45°  and 
O°<0<9O°.  The  pure-mode  eigenvalues  for  U2B4O7  are  listed  in  Table  3-1.  The  pure-mode 

/  r 

eigenvalues  (c66)  and  (c^)  correspond  to  shear  waves  propagating  through  the  plate  in  the  X2 

f 

direction.  The  pure-mode  eigenvalue  (c22 )  corresponds  to  longitudinal  waves  propagating  in  the 
x2  direction. 

The  eigenvalue  equation  (3.22)  may  now  be  solved  at  the  pure  mode  orientations  to  yield 

/  / 

eigenvalue  expressions  that  define  the  relationships  between  the  eigenvalues  (c66)  ,  (c22)  ,  and 
/ 

(c^)  and  the  material  constants.  The  eigenvalue  expressions  for  U2B4O7  are  listed  in  Tables  3- 
2  and  3-3.  Measuring  the  antiresonant  frequencies  of  crystals  cut  at  the  pure-mode  orientations 


3-10 


yields  a  set  of  pure-mode  eigenvalues  through  Equation  (3.21).  Along  with  the  eigenvalue 
expressions  of  Table  3-2  and  Table  3-3,  this  set  of  pure-mode  eigenvalues  yields  an 
overdetermined  system  of  equations  for  the  extraction  of  the  material  constants. 

Table  3-1.  Pure-Mode  Loci  in  Dilithium  Tetraborate  [2]. 


Locus  Distinct  Eigenvalues 


0°«t><45°,  e=o° 

/ 

(^44 ) 

0°<(j)<45o,  0=90° 

ft.)' 

/ 

(^22) 

ft.) 

(|)=0o,  0=0° 

(ft.) 

/ 

(^22) 

/ 

(^44) 

<h=0°,  O°<0<9O°,  0*22.7° 

ft.) 

(J)=0o,  0=22.7° 

ft*)' 

/ 

CC22  ) 

/ 

(^44) 

(|)=45°,  0=0° 

ft.)' 

/ 

(C22  ) 

60 

<J>=45°,  O°<0<9O°,  09*32.4° 

fe) 

(|)=45°,  0=32.4° 

ft.) 

/ 

(^22) 

(W 

In  order  to  measure  the  antiresonant  frequencies  of  the  pure-mode  crystal  orientations,  it 

/ 

is  necessary  to  excite  the  proper  longitudinal  waves  (corresponding  to  (c22)  )  and  shear  waves 

/  / 

(corresponding  to  (c66)  and  (c^)  ).  The  acoustic  waves  are  excited  with  electric  fields  directed 
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Table  3-2.  Symmetry  Class  Pure-Mode  Eigenvalue  Expressions  for  Dilithium 
_ Tetraborate  [2]. _ 


Locus 

Eigenvalue  Expressions 

0°«J><45°,  0=0° 

(?„)'=  4 

bll 

0°<<|)<45o,  0=90° 

(c66)  =  4 

(c  )=cE+i- 

VC22 )  C33  T  s 

E33 

(C44)  =4 

<j)=0°,  0=0° 

£ 

II 

(c22)  -c„ 

(c  1  =cE  +^- 

VC44 )  C44  T  5 

en 

<t>=0°,  O°<0<9O° 

(c66)  =cos20(4)  +  sin20(c^) 

<))=450,  0=0° 

(C6fi)  =  *2C11  ”2*  C12 

(^22)  =  2*^11  ^"'2C12  ~*"C66 

(c  )  =CE  +  flL 
\  44  /  l«T  s 

fcn 

<|)=450,  O°<0<9O° 

(c66)  =  cos2 0(^4  -2  42)+ sin2 0(4) 

through  the  thickness  of  the  crystal  plate  (thickness  excitation,  designated  TE)  or  with  electric 
fields  directed  along  the  surface  of  the  crystal  plate  (lateral  excitation,  designated  LE).  The  lateral 
field  excitation  field  may  be  applied  in  any  direction  along  the  face  of  the  crystal,  however,  in 
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Table  3-3.  Material-Specific  Pure-Mode  Eigenvalue  Expressions  for  Dilithium 
_ Tetraborate  [2], _ 


Locus 

Eigenvalue  Expressions 

<f)=0°,  0=22.7° 

(««)’=  0.85(4)+ 0.15(4) 

(c22 )  =  0.72(4  )  +  0.25(4  +  2c£  )  +  0.02(4  ) 

+[0.33(2e15  +  e3i )+  0.06(e33 )] 2  /[o.85(e  ?, )  +  0.1 5(e  33 )] 

(^44  )  =  0.49^44 )  +  0.1 3(cl !  +  c33  —  2  c13 ) 

+[0.79(e,s)+0.14(ejj -eu  -«„)]' !/[0.85(ef,)+0.15(e’,)] 

(J)=45°,  0=32.4° 

(c66 )'  =  0.36(4  )  -  0-36(4, )  +  0.29(4 ) 

(c,!)'=0.25(c'+4+24)  +  0.41(c*+2c*)+0.08(4) 

+[0.38(2els+c„)+0.15(e„)]7[0.7l(ef,)+0.29(e;,)] 

(4)  - o.i8(4)+o.io(4  +4+24+24  -44) 

+[0.60(e]!)+0.24(<.J3  -e,,)] 1  /[o.7  l(e  f, ) + 0.29(e  )] 

order  to  completely  characterize  the  crystal’s  response  it  is  only  necessary  to  define  two 
orthogonal  excitation  directions  along  the  face  of  the  crystal.  It  is  necessary  to  define  another 
rotation  angle,  \j/,  resulting  in  the  IEEE  triply  rotated  plate  designation  (YXwlt)(|>,0,V|/  [17].  The 
“t”  denotes  a  third  right-handed  rotation  about  the  doubly  rotated  thickness  axis  (x  2-axis)  of  \|/ 
degrees.  The  orientation  (YXwlt)<j>,0,O  corresponds  to  lateral  field  excitation  in  the  \j/=0° 
direction  of  a  (YXwl)<|),0  doubly  rotated  crystal.  This  orientation/excitation  combination  has  the 
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designation  LEO.  The  orientation  (YXwlt)<|),0,9O  corresponds  to  a  (YXwl)<|),0  doubly  rotated 

crystal  with  lateral  field  excitation  in  the  \|/=90°  direction.  This  orientation/excitation  combination 

has  the  designation  LE90.  Therefore,  the  only  remaining  question  is  which  modes  may  be  excited 

electrically.  This  question  is  addressed  by  Kosinski  who  concludes  that  all  but  two  of  the  pure- 

modes  may  be  excited  electrically  [2].  The  pure-modes  that  cannot  be  excited  electrically  are  the 

/ 

modes  associated  with  (c66)  for  the  (YXwl)0,0  (Y-cut)  and  (YXwl)45,0  orientations  [2].  In 

general,  modes  may  be  excited  by  both  thickness  and  lateral  excitation,  however,  some  modes 
may  only  be  excited  by  either  TE,  LEO,  or  LE90  [2].  A  sample  set  suggested  by  Kosinski  for  the 
determination  of  all  of  the  material  constants  of  U2B4O7  is  given  in  Table  3-4  [2].  Seventeen 

different  crystal  orientation/excitation  combinations  are  utilized  in  the  sample  set  [2]. 
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Table  3-4.  Suggested  Sample  Set  Summary  for  Dilithium  Tetraborate  [2]. 
Shown  are  the  orientation/excitation  combinations  recommended  for  the  extraction 
of  the  various  constants  or  combinations  of  constants. 


Constant(s) 

Orientation 

Excitation 

£S 

cll 

(YXwl)0,0 

(Y-cut) 

TE 

£S 

c33 

(YXwl)0,90 

(Z-cut) 

TE 

(YXwl)0,45 

TE 

£S  £S 
fcll  »  fc33 

(YXwl)45,45 

TE 

(YXwl)0,22.7 

TE 

(YXwl)45,32.4 

TE 

e\5 

(YXwl)0,0 

(Y-cut) 

TE 

(YXwl)45,0 

TE 

^33 

(YXwl)0,90 

(Z-cut) 

TE 

e\5  ’  ^31  *  ^33 

(YXwl)0,22.7 

TE 

(YXwlt)0, 22.7,90 

LE90 

(YXwl)45,32.4 

TE 

(YXwlt)45,32.4,90 

LE90 

4 

(YXwlt)0,0,90 

(Y-cut) 

LE90 

4 

(YXwl)0,90 

(Z-cut) 

TE 

4 

(YXwl)0,0 

(Y-cut) 

TE 

(YXwlt)0,90,0 

(Z-cut) 

LEO 

(YXwlt)0,90,90 

(Z-cut) 

LE90 

(YXwl)45,0 

TE 

E  E 

^44  *  ^66 

(YXwlt)0,45,0 

LEO 

(YXwlt)0, 22.7,0 

LEO 

E  E  E 

C||,  CI2, 

(YXwlt)45,45,0 

LEO 

(YXwlt)45,32.4,0 

LEO 

Cll»  C12>  C66 

(YXwlt)45,0,90 

LE90 

E  E  E  E 

C 11  »  ^13  >  ^33  >  C44 

(YXwl)0,22.7 

TE 

(YXwlt)0,22.7,90 

LE90 

E  E  E  E  E  E 

>  C\2  >  C13  »  C33  »  C44  ’  C66 

(YXwl)45,32.4 

TE 

(YXwlt)45, 32.4,90 

LE90 
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Chapter  4  -  Experimental  Procedure 


This  chapter  presents  the  experimental  approach  utilized  in  this  thesis.  The  first 
section  provides  an  introduction  to  the  piezoelectric  material  constants  measurement 
technique  used  in  this  thesis.  The  second  section  describes  the  equipment  and  the  crystal 
samples  as  well  as  the  experimental  setup.  The  third  section  discusses  the  experimental 
method.  The  fourth  section  describes  the  experimental  procedure. 

4.1  Introduction 

The  purpose  of  the  experimental  technique  utilized  in  this  thesis  is  to  determine  the 
values  for  the  stiffness  eigenvalues  of  multiple  samples  of  U2B4O7  over  a  range  of 

temperatures.  In  order  to  accomplish  this  goal,  a  sample  set  composed  of  U2B4O7 

resonators  cut  at  selected  pure-mode  orientations  is  utilized.  The  orientations  of  the 
sample  set  crystals  are  chosen  to  provide  an  overdetermined  set  of  eigenvalue  equations 
(see  Tables  3-2  and  3-3)  for  use  in  the  extraction  of  the  selected  material  constants  of 
U2B4O7.  Using  eigenvalue  expressions  along  with  the  determined  values  of  stiffness 

eigenvalues  for  each  selected  crystal  orientation,  the  selected  material  constants  may  be 
extracted  using  a  data  analysis  process  discussed  in  Chapter  5. 

Choosing  the  crystal  orientations  that  will  be  included  in  the  sample  set  involves 
the  analysis  of  the  eigenvalue  expressions  of  Tables  3-2  and  3-3.  Assuming  the  value  ef, 
is  known  (obtained  from  the  measurements  of  Shiosaki  [10]),  the  squared  value  of  el5  may 
be  considered  a  linear  term  in  the  symmetry  class  pure-mode  eigenvalue  equations.  Nine 
of  the  eleven  material  constants  of  U2B4O7  appear  in  the  linear  symmetry  class  pure- 
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mode  eigenvalue  equations.  The  remaining  two,  cf3  and  e3h  appear  only  in  the  nonlinear 

material-specific  eigenvalue  equations.  Three  of  the  nine  that  appear  in  linear  equations, 

/ 

4,  e33,  £33 ,  appear  in  only  one  linear  equation.  The  eigenvalue  equation  for  (c22)  for 
(YXwl)O°<<|)<45o,0=9Oo  is  the  only  linear  equation  relating  these  three  constants, 
therefore,  they  cannot  be  obtained  from  a  linear  least  squares  extraction.  The  five 
remaining  material  constants  of  U2B4O7,  c£,  cf2,  c£,  c6£6,  and  e15 ,  occur  in  multiple 

linear  eigenvalue  equations  which  form  an  overdetermined  set  of  linear  equations.  These 
are  the  selected  material  constants  whose  temperature  coefficients  will  be  determined  with 
the  use  of  a  linear  least  squares  matrix  process  derived  in  Chapter  5.  All  but  one  of  the 
symmetry  class  eigenvalue  expressions  of  Table  3-2  are  linear  equations  containing  the 
selected  material  constants.  Therefore,  all  but  one  symmetry  class  pure-mode  crystal 
orientation/excitation  combinations  may  be  used  for  the  determination  of  the  selected 
material  constants.  The  experimental  procedure  in  this  thesis  utilizes  seven  different 
crystal  orientations;  with  two  of  the  selected  orientations  excited  using  both  TE  and  LE. 
Table  4-1  shows  the  sample  set  chosen  for  this  thesis  as  well  as  the  corresponding 
eigenvalue  expressions.  When  these  orientation/excitations  are  used  to  measure  the 
stiffness  eigenvalues,  the  sample  set  gives  nine  equations  for  the  determination  of  five 
unknowns.  The  accuracy  of  the  calculated  material  constants  is  improved  because  the 
sample  set  results  in  an  overdetermined  system  of  equations.  Since  the  value  ef,  is 
obtained  from  previously  published  work,  the  value  of  ef5  may  be  obtained  through  the 
linear  least  squares  method  of  Chapter  5. 
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Table  4-1.  Sample  Set  and  Eigenvalue  Expressions. 


Number  of 

Crystal 

Samples 

Orientation  and  Excitation 

Eigenvalue  Expression 

4 

(YXwl)0,0  TE 

/  2 
(c  )  =cE 

VC44;  C44  ^  S 

tll 

4 

(YXwlt)0,0  LE90 

(c22)  =Cu 

4 

(YXwlt)0,45,0  LEO 

(%e)  =  2(c6£6)+2(c^) 

2 

(YXwl)45,0  TE 

fcd  =4+4- 
E11 

4 

(YXwlt)45,0,90  LE90 

(C22)  =2'Cll  2‘C12  '*_C66 

4 

(YXwlt)45,45,0  LEO 

)  =  j{lCU  ~2Cn)+^{Cu) 

4 

(YXwlt)0,90,90  LE90 

(%6) 

4 

(YXwlt)0,28.2,0  LEO 

(c66)'=0.78(c6£6)+0.22(4) 

4 

(YXwlt)45,56.1,0  LEO 

(c66)'=0.3l(ic£-|c1£2)  +  0.69(4) 

In  order  to  determine  the  stiffness  eigenvalues  of  the  pure-mode  crystal 
orientations,  it  is  necessary  to  measure  the  zero  mass  loading  fundamental  antiresonance 
frequency  of  the  crystal.  Knowing  mass  density,  and  thickness  of  the  crystal 
resonator  sample,  the  value  for  its  stiffness  eigenvalue  may  be  determined  through 
Equation  (3.21).  Measurements  of  the  mass  density  of  the  crystal  samples  ranged  from 
2414.6  kg/m3  to  2453.2  kg/m3  with  a  relative  accuracy  of  0.1%  [2].  Measurements  of  the 
thicknesses  of  the  samples  ranged  from  196.0  (im  to  206.9|Xm  with  a  relative  accuracy  of 
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0.1%  [2].  The  f™'s  are  measured  using  the  experimental  setup  described  in  the 
following  section.  A  complete  description  of  all  crystal  samples  and  fj$  data  is  presented 
in  Appendix  A. 

4.2  Experimental  Setup 

The  measurements  of  the  /jo’s  of  the  crystal  samples  were  performed  using  a 

Hewlett-Packard  HP  4195A  Network/Spectrum  Analyzer  and  Measurement  Unit,  an 
Anzac  Mod  HI  Hybrid  Junction,  a  General  Radio  Company  874-VCL  Variable  Capacitor, 
and  a  crystal  test  fixture.  The  crystal  test  fixture  contained  nine  Augat  8000- AG  1  crystal 
mounts  that  were  fastened  inside  a  0.25  inch  thick  aluminum  box  shown  in  Figure  4-1. 


Temperature  Probe  Ports 


H -  5”  - H 


Fig.  4-1.  Crystal  test  fixture. 
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Temperature  was  controlled  by  a  Lab-Line  Ultra-Clean  100  Oven  Model  3490M  and 
measured  with  an  Omega  Model  DP95  RTD  Thermometer  equipped  with  two  Omega 
PRP-1  High  Precision  RTD  Probes.  The  experimental  setup  is  shown  in  Figure  4-2. 


Fig  4-2.  Experimental  Setup. 

Accuracy  in  the  temperature  measurements  is  a  critical  parameter  for  extracting 
accurate  results  for  the  temperature  coefficients  of  the  material  constants.  The  issue  of 
temperature  accuracy  was  addressed  by  the  construction  of  the  test  fixture  and  the  use  of 
dual  thermometry  for  temperature  measurement.  The  test  fixture  serves  the  dual  purpose 
of  holding  the  crystals  and  stabilizing  the  temperature.  By  completely  enclosing  the 
crystals  in  a  box  of  aluminum,  the  test  fixture  stabilizes  the  temperature  of  the  crystal 
samples  in  two  ways.  First,  the  aluminum  acts  like  a  “heat  capacitor,”  storing  the  heat 
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energy  of  the  oven.  Once  the  aluminum  temperature  has  stabilized  at  the  average  oven 
temperature,  small  changes  or  instabilities  in  the  oven  temperature  have  little  effect  on  the 
aluminum  temperature  due  to  the  metal’s  high  thermal  capacity.  The  aluminum  box, 
therefore,  serves  to  stabilize  temperature  by  surrounding  its  contents  with  a  layer  of  stable 
heat  energy.  The  test  fixture  also  stabilizes  the  temperature  of  the  crystal  samples  by 
reducing  the  effect  of  convection.  By  completely  enclosing  the  crystals,  the  test  fixture 
prevents  the  thermal  convection  currents  of  the  oven  from  reaching  the  crystals. 

In  order  to  accurately  measure  the  temperature  of  the  crystals,  two  temperature 
probes  were  used.  In  addition  to  providing  two  temperature  measurements  at  each  data 
point,  this  dual  thermometry  provided  a  rough  measurement  of  thermal  gradients  in  the 
test  fixture.  The  difference  in  the  temperature  measurements  of  each  probe  was  used  in 
determining  the  stability  of  the  temperature  for  each  data  point.  Therefore,  the  use  of  dual 
thermometry  provided  more  accurate  temperature  measurements. 

4.3  Experimental  Method 

In  order  to  determine  the  stiffness  eigenvalues  of  Equation  (3.21),  it  is  necessary  to 
measure  mass  density,  thickness,  and  fundamental  zero  mass-loading  antiresonant 

frequency  (fJo)  of  the  crystal.  The  accuracy  of  the  fjg  measurements  is  another 

parameter  critical  to  the  accurate  determination  of  the  temperature  coefficients  of  the 
material  constants.  The  experimental  procedure  utilized  in  this  thesis  enhances  the 
accuracy  of  the  measurements.  The  accuracy  of  the  antiresonant  frequencies  is 

enhanced  through  the  measurement  of  multiple  harmonics.  Using  the  exact  transmission 
line  analog  for  an  infinite  piezoelectric  plate  resonator  [18-25],  it  can  be  shown  that 
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measuring  multiple  harmonics  of  the  mass  loaded  resonant  frequency  )  eliminates 


the  effect  of  piezoelectric  coupling.  For  thickness  excitation,  the  resonance  frequencies 
are  related  through  the  following  exact  transmission  line  equation  [26]: 

X 


tan(X)  = 


,  X  = 


_  fW) 
ft  J R\ 


r  Rii 


(4.1) 


*2+HX!’“  2  /« 

The  value  p  is  a  unitless  parameter  called  the  mass  loading,  is  the  zero  mass  loading 
fundamental  antiresonant  frequency,  is  the  Mth  order  mass  loaded  resonant 

frequency,  and  k  is  the  piezoelectric  coupling  constant  defined  as  follows  [2]: 


k2  = 


8  Sc 


(4.2) 


As  the  harmonic  number  increases  in  Equation  (4.1),  f^/M  converges  to  the 
fundamental  mass  loaded  antiresonant  frequency.  In  other  words,  at  higher  harmonics 
/jjJ0  becomes  an  odd  integer  multiple  of  the  fundamental  mass  loaded  antiresonant 


frequency.  This  is  shown  graphically  in  Figure  4-3. 

In  order  to  obtain  the  fundamental  zero  mass  loading  antiresonant  frequencies  to 
be  used  with  Equation  (3.21),  it  is  necessary  to  consider  the  effects  of  mass  loading.  The 
mass  loading  of  a  crystal  resonator  is  defined  as  the  ratio  of  the  mass  of  the  plated 
electrodes  to  the  mass  of  the  underlying  crystal  [27].  Mass  loading  has  the  effect  of 
lowering  the  antiresonant  frequencies  of  a  crystal  resonator.  Since  neither  the  mass  of  the 
electrodes  or  the  mass  of  the  crystal  are  functions  of  temperature,  mass  loading  is  not  a 
function  of  temperature.  Therefore,  in  order  to  eliminate  the  effects  of  mass  loading  on 
frequency  versus  temperature  curves  it  is  sufficient  to  normalize  the  curve  to  its  room 
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Figure  4-3.  Graphical  representation  of  the  transcendental  Equation  (4.1) 
with  k  =  .25  and  |i  =  .0015.  As  the  order  of  the  harmonic  increases,  the 

solutions  (circled)  approach  odd  integer  values,  thus  the  value 

approaches  M. 

temperature  zero  mass  loading  antiresonant  frequency.  This  technique  is  used  in  Chapter 
5  to  build  the  antiresonant  frequency  versus  temperature  curves  used  in  the  data  analysis. 
The  higher  order  mass  loaded  resonant  frequencies  of  the  LE  modes  converge  even  more 
quickly  to  their  corresponding  fundamental  mass  loaded  resonant  frequencies,  so 
measuring  the  s  of  LE  crystals  and  normalizing  is  also  sufficient  to  obtain  frequency 
values  for  use  in  Equation  (3.21)  [26]. 
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4.4  Experimental  Procedure 


The  experimental  procedure  involved  measurements  of  multiple  harmonics  of  the 
/W’s  of  each  crystal  sample  over  a  temperature  range  of  20°C  to  150°C.  The 

temperature  was  incremented  in  10°C  steps.  The  values  of  f^/M  converged  into  the 

fS’s.  A  measurement  trial  constitutes  a  complete  set  of  f^'s  measured  at  each 

temperature  point  for  one  crystal  sample.  There  were  53  measurement  trials  completed; 
each  crystal  in  the  sample  set  of  Table  4-1  is  measured  at  least  once.  A  complete  list  of  all 
measurements  is  given  in  Appendix  A.  Harmonics  up  to  and  including  the  ninth  harmonic 
were  measured  on  most  of  the  crystal  samples.  On  some  of  the  crystal  samples,  however, 
only  one  mode  was  measured  because  the  other  modes  could  not  be  resolved.  The  TE 
resonance  peaks  had  full  width  half  maximum  (FWHM)  values  of  approximately  5  kHz, 
while  the  LE  resonance  peaks  had  a  FWHM  of  approximately  1  kHz.  Figures  4-4  and  4-5 
show  example  TE  and  LE  resonance  peak  measurements. 

Temperature  was  measured  to  within  ±0.1  °C,  and  a  temperature  gradient 
(difference  measured  between  probes)  greater  than  the  accuracy  of  the  thermometer 

A 

occurred  only  for  the  111.2°C  (0.11°C  disparity)  and  120.9°C  (0.1 2°C  disparity) 
measurements  of  crystals  #5-12  and  #20.  With  the  fj^’s  (from  the  convergence  of 
f^/M)  measured  at  each  data  point,  the  corresponding  stiffness  eigenvalues  may  be 
calculated  from  Equation  (3.21).  The  next  step  is  data  analysis  described  in  Chapter  5. 
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BTM  □ I V  CENTER  11  672  500.  000  Hz 

-50.  00m  36.  □□  SPAN  500  □□□.  000  Hz 


Fig.  4-4.  Example  of  TE  fundamental  resonance  peak  measurement  for 
crystal  #4,  (YXwl)0,0  TE. 


DIV  DIV  CENTER  11  451  125.000  Hz 

500.  OjJ  36.  00  SPAN'  50  000.  000  Hz 


Fig.  4-5.  Example  of  LE  fundamental  resonance  peak  measurement  for 
crystal  #12,  (YXwlt)0,45,0  LEO. 
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Chapter  5  -  Data  Analysis 


This  chapter  presents  the  method  used  to  extract  the  temperature  coefficients  of 
the  selected  material  constants  of  dilithium  tetraborate.  The  first  section  is  an  overview  of 
the  data  analysis  procedure.  The  second  section  describes  the  purpose  of  the  third  order 
power  series  fit  to  the  measured  antiresonant  frequencies.  The  third  section  details  the 
calculation  of  the  stiffness  eigenvalues  from  the  antiresonant  frequency  curves.  The  third 
section  also  includes  a  comparison  of  the  calculated  room  temperature  stiffness 
eigenvalues  with  the  results  obtained  by  Kosinski  [2]  using  the  same  crystals.  The  fourth 
section  describes  the  linear  least  squares  extraction  process  for  the  material  constants. 
The  fifth  section  describes  the  calculation  of  the  temperature  coefficients  of  the  material 
constants.  The  sixth  section  is  a  summary  of  the  results. 

5.1  Data  Analysis  Overview 

Now  that  the  mass  loaded  antiresonant  frequencies  have  been  obtained  for  each 
crystal  orientation  using  the  experimental  procedure  of  Chapter  4,  the  first  step  in  the  data 
analysis  is  the  normalization  of  the  frequency  versus  temperature  curves.  The  next  step  is 
the  calculation  of  the  stiffness  eigenvalues  using  Equation  (3.21).  Both  the  mass  density 
and  the  thickness  of  the  crystals  are  functions  of  temperature,  so  each  must  be  calculated 
over  the  appropriate  temperature  range.  The  resulting  room  temperature  stiffness 
eigenvalues  were  compared  with  those  of  Kosinski  [2]  in  order  to  check  the  data 
extraction  program.  The  next  step  is  a  linear  least  squares  fit  of  the  stiffness  eigenvalue 
data.  Using  the  eigenvalue  expressions  found  in  Table  3-2  for  the  selected  sample  set 
orientations  results  in  an  overdetermined  system  of  equations  relating  stiffness  eigenvalue 
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data  to  the  material  constants.  The  linear  least  squares  process  is  used  for  the  extraction 
of  the  material  constants  from  the  overdetermined  system.  The  accuracy  of  the  values  for 
the  material  constants  is  enhanced  by  using  an  overdetermined  system  of  equations  in  the 
extraction  process.  The  data  extraction  was  performed  using  MathCAD  Plus  5.0  [28]. 
The  data  extraction  program  is  given  in  Appendix  A. 

5.2  Third  Order  Power  Series  Fit  for  Antiresonant  Frequencies 

Once  the  antiresonant  frequencies  were  measured,  it  was  necessary  to  construct 
antiresonant  frequency  versus  temperature  curves  so  that  the  data  could  be  used  in  the 
extraction  process.  A  third  order  power  series  with  a  reference  temperature  of  25°C  (see 
Equation  (2.6))  was  used  to  fit  the  antiresonant  frequency  data  for  each  trial.  Because  the 
crystal  samples  had  electrodes  plated  to  their  surface,  the  antiresonant  frequency  data  was 
normalized  to  the  room  temperature,  unelectroded  measurements  of  Kosinski  [2].  The 
maximum  deviations  of  the  measured  fundemental  mass  loaded  antiresonant  frequencies 
from  Kosinski’s  values  were  only  1.3%  for  TE  measurements  and  0.6%  for  LE 

measurements.  The  unelectroded  measurements  of  Kosinski  were  used  to  ensure 

zero-mass  loading  conditions.  Replacing  the  zeroth  order  fundemental  mass  loaded 
antiresonant  frequency  temperature  coefficients  with  the  unelectroded  measurements  of 
Kosinski  should  completely  eliminate  mass  loading  effects  in  this  thesis’  fundemental  mass 
loaded  antiresonant  frequency  data.  The  third  order  power  series  fit  provided  an  easy 
opportunity  to  normalize  the  measured  antiresonant  frequencies.  Once  the  temperature 
coefficients  of  antiresonant  frequency  were  calculated,  the  zeroth  order  temperature 
coefficient  was  replaced  with  the  room  temperature,  unelectroded  measurements  of 
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Kosinski.  The  resulting  normalized  antiresonant  frequency  versus  temperature  curves 
were  used  in  the  calculation  of  the  stiffness  eigenvalues  for  each  crystal  orientation  at  the 
temperature  points:  20°C,  30°C,  40°C,...  150°C. 

5.3  Stiffness  Eigenvalue  Determination 

Once  the  antiresonance  frequency  versus  temperature  curves  have  been 
determined,  the  next  step  is  the  determination  of  the  stiffness  eigenvalues  for  each  crystal 
orientation.  The  stiffened  elastic  constants  are  calculated  at  each  temperature  point. 

In  order  to  calculate  the  stiffness  eigenvalues  from  the  antiresonant  frequency  data, 
Equation  (3.21)  is  used.  Equation  (3.21)  shows  that  the  stiffness  eigenvalues  are 
functions  of  zero  mass  loading  fundamental  antiresonant  frequency,  mass  density,  and 
crystal  thickness.  Both  the  mass  density  and  the  thickness  of  the  crystal  samples  change 
with  temperature,  therefore,  it  is  necessary  to  calculate  their  values  at  the  appropriate 
temperature  points. 

In  order  to  calculate  the  values  of  mass  density  and  crystal  thickness  at  the 
temperature  points,  previously  reported  values  of  the  coefficients  of  thermal  expansion  of 
U2B4O7  were  used.  The  first  and  second  order  coefficients  of  thermal  expansion  were 

obtained  from  Shiosaki,  et.  al.  [10]  and  are  listed  in  Table  5-1.  Shiosaki,  et.  al.,  uses  a 
reference  temperature  of  25°C  for  the  coefficients  of  thermal  expansion  of  U2B4O7  [10]. 

The  thermal  expansion  coefficients  for  the  unrotated  Xi  direction  ((Xu’s)  are  identical  to 
the  thermal  expansion  coefficients  for  the  unrotated  x2  direction  (a22’s).  The  first  and 


5-3 


Table  5-1.  First  and  Second  Order  Coefficients  of  Thermal  Expansion  of 
U2B4O7  Obtained  from  Shiosaki,  et  al.  [10]. 


First  Order  (xlO'6/K) 

Second  Order  (xlO'9/K2) 

au 

11.1 

5.6 

&33 

-3.74 

21 

second  order  temperature  coefficients  of  mass  density  were  calculated  with  Equation 
(5.1). 

7]p<")=-«)+a£)+ajS))  (5.1) 

The  value  7’p(n)  is  the  nth  order  temperature  coefficient  of  mass  density  and  a£°  are  the 

nth  order  thermal  expansion  coefficients.  The  mass  density  of  the  crystal  samples  was 
calculated  using  a  second  order  power  series.  Equation  (5.2). 

p  =  p 0 [l  +  Tp (1) (6  -  25)  +  Tp (2) (0  -  25 )2  ]  (5.2) 

The  value  p  is  the  mass  density,  p0  is  the  room  temperature  (25°C)  mass  density,  and  0 

is  the  temperature.  The  room  temperature  values  of  mass  density  for  the  crystal  samples 
were  measured  by  Kosinski  to  an  accuracy  of  0.1%  [2]. 

Because  the  coefficients  of  thermal  expansion  of  IJ2B4O7  depend  on  crystal  cut 

orientation,  calculating  the  thicknesses  of  the  crystal  samples  at  each  temperature  point  is 
more  complicated.  In  order  to  use  the  coefficients  of  thermal  expansion  of  Table  5-1,  they 
must  first  be  transformed  to  the  appropriate  crystal  orientation  using  Equation  (5.3)  [16]. 

[a,,<n)  ]  =  [a]  [a (n)]  [a] T  (5.3) 
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The  matrix  [a]  is  the  direction  cosine  matrix  defined  in  Chapter  3,  and  matrix  [a(n)]  is  a 
matrix  of  the  coefficients  of  thermal  expansion  defined  in  Equation  5.4. 


[oc(n)]  = 


«M  0 

a(n) 
uu 


0 

0 


0 

0 


0  a 


(«) 

33 


(5.4) 


When  [a(n)]  is  transformed  using  Equation  (5.3)  for  a  given  crystal  cut  orientation,  the 
resulting  matrix  may  have  off-diagonal  elements  as  shown  in  Equation  (5.5). 


[a,,(n)]= 


Vy  "(rt) 
Uu 

f y  "( n ) 
U12 

yy/'(n) 

UI3 

n  "(n) 
u'21 

n  "(") 
^22 

n  "(«) 
IX  23 

M  "(«) 

**31 

IX  32 

yy  "(n) 

^33 

(5.5) 


Since  the  crystal  samples’  thicknesses  are  directed  along  the  x"  axis,  the  element  is 
the  nth  order  coefficient  of  thermal  expansion  for  the  thickness  of  the  doubly  rotated 
crystal  plate.  Using  the  value  fora"2(n) ,  the  thickness  of  a  crystal  sample  at  the 
appropriate  temperature  points  is  calculated  with  a  second  order  power  series.  Equation 
(5.6). 

/  =  /o[l+a''2(,)(0  -25)+a"(2)(6  -25)2]  (5.6) 

In  Equation  (5.6),  l  is  the  crystal  thickness  and  l0  is  the  room  temperature  (25°C)  crystal 

thickness.  The  room  temperature  crystal  thicknesses  of  the  investigated  crystal  samples 
were  measured  by  Kosinski  to  within  0.1%  [2]. 

Once  the  values  of  mass  density  and  crystal  thickness  are  calculated  over  the  range 
of  temperatures,  the  next  step  in  the  data  analysis  is  to  use  Equation  (3.21)  to  calculate  the 
stiffness  eigenvalues  at  the  temperature  points  for  each  measurement  trial.  Table  5-2 
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shows  a  comparison  between  the  calculated  room  temperature  values  for  the  stiffness 
eigenvalues  and  those  obtained  by  Kosinski  using  the  same  crystal  samples  [2], 

Table  5-2.  A  Comparison  of  the  Calculated  Values  of  Room  Temperature 
c  with  the  Values  Obtained  by  Kosinski  Using  the  Same  Crystals  [2]. 


Crystal  Orientation/ 
Excitation 

Results  of 
Calculated  c  ’s 
(xlO9  Pa) 

this  thesis: 

Standard  Deviation 
(%) 

Kosinski  c  ’s 
(xlO9  Pa) 

Error 

(%) 

(YXwl)0,0  TE 

59.21 

0.29 

59.12 

0.155 

(YXwlt)0,0,90  LE90 

135.95 

0.23 

135.9 

0.038 

(YXwlt)0,45,0  LEO 

52.46 

0.23 

52.49 

0.067 

(YXwl)45,0  TE 

59.00 

0.29 

58.95 

0.079 

(YXwlt)45,0,90  LE90 

115.54 

0.28 

115.4 

0.117 

(YXwlt)45,45,0  LEO 

62.61 

0.30 

62.49 

0.194 

(YXwlt)0,90,90  LE90 

57.04 

0.40 

57.06 

0.031 

(YXwlt)0, 28.2,0  LEO 

49.757 

0.27 

49.76 

0.006 

(YXwlt)45,56.1,0  LEO 

60.48 

0.25 

60.49 

0.016 

The  utility  of  Table  5-2  is  merely  to  check  the  validity  of  the  data  analysis  program  used  in 
this  thesis.  The  errors  in  Table  5-2  may  be  attributed  to  the  fact  that  Kosinski  used 
slightly  different  weightings  for  each  crystal  sample  used  in  calculating  the  c ’s. 

The  values  of  stiffness  eigenvalues  are  the  pure-mode  stiffened  elastic  constants 


_  '  _  '  _  'A 

C66  )  ’  if' 22  )  ’  (C44  ) 


derived  in  Chapter  3  and  listed  in  Table  4-1. 


The  eigenvalue 
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expressions  in  Table  4-1  constitute  an  overdetermined  system  relating  the  known  values  of 
the  stiffened  elastic  constants  to  the  unknown  values  of  the  material  constants 
(  e2) 

cH>ci E2’c44’c66’-f'  •  A  linear  least  squares  extraction  process  is  used  to  extract  the 

^  en  J 

values  of  the  material  constants  from  this  overdetermined  system. 

5.4  Linear  Least  Squares  Extraction  of  Material  Constants 

Every  measured  antiresonant  frequency  produces  a  calculated  stiffened  elastic 
constant  value.  Each  stiffened  elastic  constant  value  corresponds  to  one  of  the  eigenvalue 
expressions  of  Table  4-1.  In  this  thesis  nine  different  crystal  orientation/excitation 
combinations  were  used  as  shown  in  Table  4-1.  There  were  53  measurements  of 
antiresonant  frequency  made  at  each  temperature  point..  Therefore,  at  each  temperature 
point  there  are  53  equations  relating  five  unknowns.  In  order  to  extract  the  values  of  the 
five  unknowns,  a  linear  least  squares  method  is  used.  The  data  extraction  method  used  in 
this  thesis  follows  the  linear  least  squares  method  derived  by  Nye  [7]. 

(  e2  \ 

Each  stiffened  elastic  constant  is  a  function  of  five  values  c^,c^,c^,,c6£6,— ^ 

l  e.J 

and  thus  may  be  written  as  shown  in  Equation  (5.7). 

e2 

C  =  b\Cl |  +  b2Cl2  +  ^3C44  "*■  ^4C66  ^5  ~~S~  (5.7) 

EI1 

The  coefficients  (bn  ’s)  are  constants  that  depend  on  the  crystal  orientation/excitation.  For 

example,  the  eigenvalue  equation  of  the  (YXwlt)45,0,90  orientation  with  LE90  excitation 
has  coefficient  values  of  Vz,  Vi,  0,1,0  respectively.  In  general  there  may  be  N  values  of 
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c  measured  (N  =  53  for  this  thesis)  resulting  in  a  system  of  N  equations  as  shown  in 
Equation  (5.8). 


C,  —  bnCn  +  bnCn  +  ^13C44  ■*"  ^UC66  +  ^15  5 


^2  —  ^21^11  "f  ^22^12 


(5.8) 


F  =/?  c£+h  cE+h  cE +h  cE  +h  — 

^t/Af2c12  '  u NV~  44  ~  t/W4c66  T  ^JVS  s 


Each  measurement  results  in  the  true  value  of  c  plus  some  error,  v .  Rewriting  Equation 
(5.8)  to  include  the  error  results  in  the  following: 


v,  —  b\\cu  +  bnci2  +bi3c44  +bHc66  +bi5  .  c, 

en 

V2  =  ^2IC11  "*"^22C12 


—  bmci\ 


N2C12  +  ^V3C44  ^N4C66  +  ^N5 


(5.9) 


Equation  (5.9)  is  rewritten  in  matrix  form  to  obtain  Equation  (5.10). 

vi=bijZj-Cj  (5.10) 

where,  Zj  = 

Following  the  principle  of  least  squares,  the  value  of  zi  should  be  such  that  the  sum  of  the 
squares  of  the  errors  is  a  minimum.  This  condition  is  given  by  Equation  (5.12). 
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(5.12) 


|-(v,2  +  v22  +  v32  +•  •  •+ v* )  =  ±(vt  •  v(. )  =  0 


Using  the  product  rule,  Equation  (5.12)  may  be  written  as  Equation  (5.13). 


*%**>%-« 


*av° 


(5.13) 


By  taking  the  derivative  with  respect  to  Zj  of  both  sides  of  Equation  (5.10)  and  then 


multiplying  through  by  v, ,  Equation  (5.14)  is  obtained. 


8v,/  _  d(bijzi)  dcj/ 

Vi  A=V'^--V‘  A 


(5.14) 


Since  c ■  is  a  vector  of  measured  constants,  its  partial  derivative  is  equal  to  zero. 
Applying  the  result  of  Equation  (5.13)  yields  the  following: 


vAj  =  0 


(5.15) 


Now  multiplying  Equation  (5.10)  by  btJ  results  in  Equation  (5.16). 


VA= 


(5.16) 


Applying  the  result  of  Equation  (5.15)  and  solving  for  Zj  gives  Equation  (5.17). 


^-(bfaVblc, 


(5.17) 
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Equation  (5.17)  is  the  defining  equation  for  the  linear  least  squares  method  of  extracting 
the  material  constants  from  N  stiffened  eigenvalue  measurements.  Substituting  Equation 
(5.17)  into  Equation  (5.10)  gives  an  equation  for  the  accuracy  of  the  measurements: 


v.  =  (bJibik )  'b]i-  l)cj  (5-18) 

where  I  is  the  identity  matrix. 

5.5  Determination  of  Material  Constants  and  Their  Temperature  Coefficients 

The  material  constants  extraction  involves  the  construction  of  the  N  by  5  matrix  of 
coefficients,  btj .  Performing  the  appropriate  matrix  operations  on  by  and  multiplying  by  a 

vector  containing  the  corresponding  stiffened  elastic  constant  measurements  produces  zk , 

(  e2) 

a  vector  composed  of  material  constants  c/j,c^,c£,cf6,— f-  .  In  order  to  calculate  the 

l  eiJ 

value  of  e15 ,  it  is  necessary  to  know  the  value  of  ef, .  The  room  temperature  value  of  e f, 

was  obtained  from  Kosinski  [2]  and  its  corresponding  first  and  second  order  temperature 
coefficients  were  obtained  from  Shiosaki,  et.  al.  [10],  and  are  listed  in  Table  5-3. 

Table  5-3.  Values  for  Room  Temperature  Dielectric  Constant  of  U2B4O7 
(ef,  )  and  First  and  Second  Order  Temperature  Coefficients 


(TefI(1),refI(2))  [10]. 


<  (10'12F/m) 

Tef,0’  (10'6/K) 

Te?{2>  (W9/K2) 

76.69 

97.1 

2800 
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Values  for  ef,  at  the  appropriate  temperature  points  are  calculated  using  a  second  order 
power  series  equation  identical  to  Equation  (5.6). 

The  material  constants  extraction  is  repeated  at  each  of  the  14  temperature  points 
(20°C  through  150°C)  to  obtain  values  for  the  material  constants  over  the  temperature 
range.  The  values  of  ei5  are  calculated  from  the  fifth  element  of  the  zk  vectors  using  the 

ef,  values  at  the  respective  temperature  points.  In  order  to  obtain  the  temperature 
coefficients  of  the  material  constants,  a  third  order  power  series  is  fit  to  the  calculated 
material  constants  values. 

5.6  Summary  of  Experimental  Results 

The  zeroth  through  second  order  temperature  coefficients  of  the  material  constants 
of  U2B4O7  were  obtained  from  the  third  order  power  series  fit  to  the  material  constants. 

The  zeroth  order  (room  temperature  values)  serve  as  a  check  for  the  data  extraction 
program.  These  values  are  compared  to  the  room  temperature  material  constants  obtained 
by  Kosinski  [2]  and  other  researchers  in  Table  5-4.  The  results  obtained  for  the  first  and 
second  order  temperature  coefficients  of  the  selected  material  constants  of  U2B4O7  are 

listed  in  Table  5-5. 
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Table  5-4.  Comparison  of  the  Calculated  Selected  Material  Constants  with 
those  Determined  by  Kosinski  [2]  and  Others  at  Room  Temperature. 


Results 


of  this 
Thesis 

[2] 

[12] 

[29] 

[30] 

[31] 

[32] 

[33] 

4  (GPa) 

135.955 

135.823 

135.27 

126.7 

127.1 

135 

135.2 

127.6 

4  (GPa) 

-0.288 

-0.285 

0.109 

0.5 

0.6 

3.57 

0.8 

8.9 

4  (GPa) 

57.057 

57.072 

57.39 

55.0 

53.8 

58.5 

55.9 

57.5 

4  (GPa) 

47.700 

47.680 

47.38 

46.0 

57.4 

46.7 

47.3 

48.2 

eX5  (C/m2) 

.4017 

0.3918 

0.36 

0.36 

0.278 

0.472 

0.35 

0.39 

Table  5-5.  Calculated  First  and  Second  Order  Temperature  Coefficients  of 
_ Selected  Material  Constants  of  L^I^O?. 


First  Order  (xlO'6/K) 

Second  Order  (xl0'9/K) 

4 

-70.26 

-193.1 

4 

-67310 

114939 

4 

-14.53 

-398.1 

4 

-418.6 

504.4 

e\s 

-1320 

5163 

Figures  5-1  through  5-5  are  plots  of  the  calculated  material  constants  and  their 
corresponding  third  order  power  series  fits. 
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57.1 


57 

-E  56.9 
C44 

(GPa)  56.8 

56.7 

56.6 


20  40  (0  80  100  120  140  160 

Temperature  (°C) 

Fig.  5-3.  Plot  of  4  and  its  power  series  fit. 
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Fig.  5-4.  Plot  of  4  and  its  power  series  fit. 
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Chapter  6  -  Conclusions  and  Recommendations 

This  chapter  presents  a  summary  of  the  experimental  results  as  well  as 
recommendations  for  further  study.  The  first  section  presents  a  comparison  of  the 
calculated  temperature  coefficients  of  the  selected  material  constants  of  U2B4O7  to 

previously  published  values.  The  second  section  discusses  how  the  knowledge  gained  in 
this  thesis  might  be  applied  and  improved  upon. 

6.1  Comparison  of  Results 

The  first  and  second  order  temperature  coefficients  of  the  material  constants  cf, , 
4 ,  4 ,  4 ,  and  e15  of  U2B4O7  have  been  calculated  using  an  improved  resonator 
method  described  by  Kosinski  [26].  Table  6-1  compares  the  calculated  first  order 
temperature  coefficients  with  previously  published  values. 

Table  6-1.  Comparison  of  Calculated  First  Order  Temperature  Coefficients 


(x10~6/K)  with  Previously  Published  Values. 


Results  of 
this  thesis 

[29] 

[30] 

[31] 

[12] 

[13] 

4 

-71.26 

-125 

-51 

-81 

-80 

-65 

4 

-67310 

14000 

1600 

3370 

20000 

-8390 

4 

-14.53 

-23 

-22 

-18 

13 

108 

4 

-418.6 

-480 

-200 

-272 

-480 

-220 

e\5 

-1320 

-1300 

-349 

-1050 
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Table  6-2  presents  a  comparison  of  the  second  order  temperature  coefficients  with  those 
of  Shiosaki,  et.  al.  [34] 


Table  6-2.  Comparison  of  Second  Order  Temperature  Coefficients  (xlO'9/K) 
_ with  Those  of  Shiosaki,  et.  al.  [34]. _ 


Results  of 
this  thesis 

[34] 

4 

-193.1 

-440 

4 

-114939 

-17400 

4 

-398.1 

500 

4 

504.4 

-450 

e\5 

5163 

-900 

The  calculated  values  of  the  temperature  coefficients  of  the  selected  material 
constants  were  used  to  predict  the  measured  data.  This  serves  as  a  check  for  the 
calculated  results.  When  designing  piezoelectric  plate  resonators  and  transducers,  it  is 
desirable  that  the  /j J/’s  are  known  with  an  accuracy  of  at  least  1%  [2].  The  mean  error 

between  the  measured  and  calculated  /jo’s  was  1.20%  with  a  standard  deviation  of 

0.005%  for  TE  and  0.285%  with  a  standard  deviation  of  0.031%  for  LE.  The  difference 
in  the  accuracies  of  the  predictions  is  attributed  to  the  fact  that  the  measured  higher 
harmonic  resonant  frequencies  of  TE  modes  converge  more  slowly  than  the  higher 
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harmonic  resonance  frequencies  of  LE  modes.  Appendix  A  contains  the  complete 
comparison  between  all  measured  and  calculated  values  of  ■ 

6.2  Recommendations  for  Further  Study 

The  knowledge  gained  in  this  thesis  should  be  useful  in  the  development  of 
U2B4O7  based  electronics.  The  most  desirable  benefit  obtained  from  the  results  of  this 

thesis  will  probably  be  the  prediction  of  new  zero  temperature  coefficient  orientations  of 
U2B4O7.  These  new  zero  temperature  coefficient  orientations  of  Ii2®4®7  might  also 
have  advantages  over  current  zero  temperature  coefficient  cuts.  In  order  to  test  the 
validity  of  the  results  presented  here,  zero  temperature  coefficient  cuts  of  U2B4O7  may 

be  predicted  using  the  calculated  temperature  coefficients  and  tested  in  later  research. 

Further  studies  of  U2B4O7  should  focus  on  the  temperature  coefficients  for  the 

remaining  material  constants.  Also,  the  temperature  coefficients  of  expansion  and 
dielectric  constant  should  be  measured.  Accomplishing  these  two  research  goals  will 
allow  the  complete  thermal  characterization  of  U2B4O7  without  the  need  of  previously 

published  material. 

Another  possible  area  of  further  research  is  below-room-temperature 
measurements  of  U2B4O7.  The  results  of  this  thesis  are  limited  to  roughly  20°C  to 

150°C.  The  power  series  fits  that  were  used  to  obtain  the  temperature  coefficients  of  the 
selected  material  constants  should  be  most  accurate  in  the  middle  of  the  temperature  range 
and  least  accurate  near  the  extremes  of  the  temperature  range.  In  order  to  improve  the 
overall  accuracy  of  the  calculated  temperature  coefficients,  especially  at  room  temperature 
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values,  measurements  should  be  made  of  the  below-room-temperature  characteristics  of 
the  material  constants  of  026407. 
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Crystal  Sample  Number  Orientation  Excitation  Thickness  (pm)  Mass  Density  (kg/m3)  (MHz)  [35] 
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Crystal  Sample  Number  Orientation  Excitation  Thickness  (pm)  Mass  Density  (kg/m3)  fU)A0  (MHz)  [35] 
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Now  plot  the  data  and  the  power  series  fit  on  the  same  graph.  Check  the  fit  using  Chauvenet’s  criterion. 
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Now  combine  the  temperature  coeff  of  antires.  frequency  data  with  the  values  for  room  temperature  antires. 
frequency  of  Kosinski  [34]  and  calculate  the  stiffness  eigenvalues  for  each  crystal  cut  at  various  temperatures: 


Fitl8  :=  1 1.87 1  +  Coefl8  y  10'6  (j- 10-  5)  +  Coefl8 1 10_9  (j-10-  5)2  +  Coefl8  y  10~12(j- 10-  5)3 
Fit  18b.  :=  1 1.87-f  1  -t-Coefl8b  10"6  (j  10- 5)  +  Coefl8b  z  10'9(j- 10-  5)2+  Coefl8b  3  10' 12(j- 10-  5) 
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Fit40.  =12.001  jl  +  Coef40  ,10'6  (j  10-  5)  +  Coef40  2  10_9(j- 10-  5)2+ Coef403  10‘12(jl0-  5)3 
Fit49.  =  1 1.538  1  +  Coef49  j  l0'6  (j  l0-  5)+ Coef492  10'9(jl0-  5)2+ Coef49  3- 10‘12(j-10-  5)3 
Fit49b.  =  1  1.538  1  4-  Coef49b  j  106  (j  10- 5)  +  Coef49b  2-  10'9(j-10-  5)2-i-  Coef49b  y  1012-(j- 10-  5)3 


12.265  [1  +-  Coef64b  ,10'6  (j  l0-  5)+ Coef64b2-10'9  (j  l0-  5)2+ Coef64b  3- 10’12  (j-10-  5)3 


Mass  density,  crystal  thickness,  and  dielectric  constant  all  vary  with  temperature.  Therefore,  they  must  be  calculated  for  each  temperature. 
Crystal  thickness  is  also  a  function  of  orientation,  so  the  proper  rotations  must  also  be  applied. 

Mass  density  changes  with  temperature,  but  does  not  change  with  crystal  orientation: 


A-10 


For  each  of  the  30  crystals  in  the  sample  set,  the  densities  over  our  temperature  range  are: 


p21  2436.3*^  i  +■  p  1U-J  -  !>)  +  p  2*(  1U  J  -  5) 
p22.  :  =  2437.0*  1  +  p  r(10j  -  5)  +  p2  (10 j  -  5): 


Shiosaki  s  values  for  the  coefficients  of  thermal  expansion  [10]  must  be  transformed  using  the  following  rotation  matrix: 


The  first  order  coeff  is:  aO_45  ,  =  3.68  I0'6  The  second  order  coeff  is:  a0_45  2  :  =  1.32- 10 

The  thickness  as  a  function  of  temperature  is:  t0_45.  :=  202.2*  10“6  1  + a0_45  ^(lOj-  5)  +  a0_45  2-(  10  j  -  5)2] 


For  the  (YXwl)45,0  orientation: 


For  the  (YXwl)0,28.2  orientation: 
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Dielectric  constant  varies  with  temperature  and  crystal  orientation,  however,  only  the  unrotated  dielectric 
constants  are  needed,  so  no  rotations  need  to  be  performed: 
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Now  a  third  order  power  fit  is  performed  so  the  room  temperature  values  of  the  stiffened  elastic  constants  can  be  compared  to  those  of  Kosinski  [2]: 
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This  is  a  plot  of  the  cbarl  values  and  the  corresponding  power  series  fit: 


temps- 


fitcbar9  :=  genfit( temps, cbar9,  guess, Fitfunction)  -1 .944*  10 

fitcbar9  = 

-1.019*10 
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Linear  Least  Squares  Material  Constants  Extraction: 
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Our  linear  least  squares  extraction  matrix  (LLS)  is: 


TempMatrixl  - augment(cbarl,cbarlb)  TempMatrix2  :  =  augment( TempMatrixl ,cbar2)  TempMatrix3  :=  augment(TempMatrix2,cbar2b) 

TempMatrix4  -  augment(  TempMatrix3 ,  cbar3 )  T empMatrix5  :  =  augment(  T empMatrix4 ,  cbar4 )  T empMatrix6  :  -  augment(  T empMatrixS ,  cbar5 ) 

TempMatrix7  -  augment(  TempMatrix6 ,  cbar6 )  T empMatrix8  :  =  augment^  T empMatrix7 ,  cbar7 )  T empMatrix9  :  =  augment(  T empMatrix8 ,  cbar8 ) 
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Perform  the  linear  least  squares  matrix  operation: 
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The  maximum  percent  error  of  the  calculated  stiffened  eigenvalues:  MaximumError  =  0.687 
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The  Temperature  Coefficients  of  the  Material  Constants  of  Dilithium  Tetraborate: 


1.34*10 


The  standard  deviation  for  the  cEn  power  series  fit  and  corresponding  Chauvenet's  criteria: 

N  :=  14  i  :=0..N-  1 


O 


Qh 
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The  standard  deviation  for  the  cE12  power  series  fit  and  corresponding  Chauvenet's  criteria: 
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temps- 


The  standard  deviation  for  the  cE44  power  series  fit  and  corresponding  Chauvenet's  criteria: 


o 
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4.5*10 


The  standard  deviation  for  the  cEgg  power  series  fit  and  corresponding  Chauvenet's  criteria: 
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fite  15:=  genfit(  temps ,  e  1 5 ,  guess ,  F  itfunction ) 


temps-. 


The  standard  deviation  for  the  e15  power  series  fit  and  corresponding  Chauvenet's  criteria: 
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Prediction  of  Measured  Antiresonant  Frequencies  Using  Calculated 
Temperature  Coefficients  of  Selected  Material  Constants: 

Measured  Data: 
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The  following  is  a  plot  of  the  measured  frequency  curve  and  the  predicted  frequency  curve  for  Trial  #1 ,  crystal  #1 : 
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Measured  Data: 
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The  following  is  a  plot  of  the  measured  frequency  curve  and  the  predicted  frequency  curve  for  Trial  #5,  crystal  #5: 
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For  Trial  #38,  crystal  #38,  (YXwl)0,90  LE90: 
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For  Trial  #39,  crystal  #39,  (YXwl)0,90  LE90: 
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For  Trial  #40,  crystal  #40,  (YXwl)0,90  LE90: 
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Measured  Data: 
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Measured  Data: 
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